Nucleon spin densities in a light-front constituent quark model 



B. Pasquini and S. Boffi 

Dipartimento di Fisica Nucleare e Teorica, Universita di Pavia, and 
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

(Dated: February 3, 2008) 

The first moment of longitudinal and transverse spin densities of quarks in the 
nucleon are calculated in a light-front constituent quark model for the different cases 
of quark and nucleon polarization. Significant distortions are found for the transverse 
spin densities. In particular the Sivers function is predicted with opposite sign for up 
and down quarks and the Boer-Mulders function is predicted large and negative for 
both up and down quarks, in agreement with lattice calculations. Quite a different 
spin distribution is obtained for up and down quarks in the cases of quarks and 
proton transversely or longitudinally polarized in the same direction. 
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At leading order in the deep inelastic regime the quark structure of the nucleon is fully 
parameterized in terms of generalized parton distributions (GPDs) (for recent reviews, see 
Refs. P, El E])- Four of them, usually indicated by H, E, H, E, are helicity conserving 
(chiral even) and the other four, denoted by Ht, Ht, E^, Et, are helicity flip (chiral 
odd) distributions [H E]. They are all functions of three kinematical variables, i.e. the 
longitudinal momentum fraction x of the quark, the invariant momentum transfer t and 
the skewness parameter ^ describing the fraction of the longitudinal momentum transfer to 
the nucleon. For each flavour they are related to important quantities characterizing the 
nucleon structure. In the forward limit, i.e. t = 0, H, H and reduce to the parton 
density fi{x), helicity distribution gi{x) and transversity distribution hi{x), respectively. 
First moments of H, E, H, and E are identified with the quark contribution to the Dirac 
and Pauli form factors, Ei and E2, and the axial and pseudoscalar form factors, Ga and Gp, 
respectively. Thus in the forward limit E and H are related to the quark contribution to the 
anomalous magnetic moment k and the axial- vector coupling g^. The first moment of the 
combination Et + 2Ht in the forward limit gives the quark contribution to the anomalous 
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tensor magnetic moment kt- Finally, Et vanishes for ^ = 0. 

When = and x > 0, by a two-dimensional Fourier transform to impact parameter 
space GPDs can be interpreted as densities of quarks with longitudinal momentum fraction 
X and transverse location h with respect to the nucleon center of momentum [HI E]. De- 
pending on the polarization of both the active quark and the parent nucleon, according to 
Refs. [TJ IE] one defines three-dimensional densities p(x, b. A, A) and p{x, b, s, S) representing 
the probability to find a quark with longitudinal momentum fraction x and transverse posi- 
tion b either with light-cone helicity A (= ±1) in the nucleon with longitudinal polarization 
A (= ±1) or with transverse spin s in the nucleon with transverse spin S. They read 



p{x,b,X,A) 



H{x, b^) + Ve^'S'^ E\x, b^) + \AH{x, b^) 



(1) 



p{x,b,s,S) 



H{x,b') + s'S' [HT{x,b- 



4M2 



AbHT{x,b' 



+ 



M 



(S'E'{x, b^) + s' [E!r{x, b^) + 2H't{x, b^ 



+ s\2bV -bH,,)S^^^H'^{xy) 



(2) 



The distributions H, E, H, Ht, etc. are the Fourier transform of the corresponding GPDs, 
i.e. 



f{x,b' 



i6A 



/(x,^ = 0,t = -A2), 



(3) 



where A is the transverse momentum transfer to the nucleon. As such, they depend on b 
only via its square b^ = 6^ thanks to rotation invariance. In Eqs. ([l| and ^ the shorthand 
notations 



_d_ 



/' / (^^2) 



d_d_ 



(4) 



have been used, and the two-dimensional antisymmetric tensor e^^ has been introduced with 



.12 



— e:^^ = 1 and e 
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= 0. M is the nucleon mass, and Roman indices are to be 



summed over. 

In Eq. ([T| the first term with H describes the density of unpolarized quarks in the 
unpolarized proton. The term with E' introduces a sideways shift in such a density when 
the proton is transversely polarized, and the term with H refiects the difference in the density 
of quarks with helicity equal or opposite to the proton helicity. 
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In the three hnes of Eq. ^ one may distinguish the three contributions corresponding 
to monopole, dipole and quadrupole structures. The unpolarized quark density in the 
monopole structure is modified by the chiral-odd terms with Ht and AbHT when both the 
quark and the proton are transversely polarized. Responsible for the dipole structure is either 
the same chiral-even contribution with E' from the transversely polarized proton appearing 
in the longitudinal spin distribution ([!]) or the chiral-odd contribution with E'j, + 2H!p from 
the transversely polarized quarks or both. The quadrupole term with is present only 
when both quark and proton are transversely polarized. 

Lattice calculations accessing the lowest two x-moments of the transverse spin densities 
of quarks in the nucleon have recently been presented [9], and impact parameter depen- 
dent parton distributions in phenomenological models of hadron light-cone wave functions 
(LCWFs) have been studied in Ref. [ID]. 
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Figure 1: The monopole contribution (left) for unpolarized quarks, the dipole contribution 
— \sxhy{E'rp + 2H'rp)/M (middle) for (transversely) x-polarized quarks, and the sum of both (right) 
in an unpolarized proton. The upper (lower) row gives the results for up (down) quarks. 



In this Letter the first x-moments of the spin distributions 

p(b, A,A)= J dx p{x,b, X,A), p{b,s,S) = J dx p{x,b, s, S) (5) 
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Figure 2: With unpolarized quarks the monopole contribution (left) for an unpolarized proton, 
the dipole contribution —^SxbyE'/M for a (transversely) x-polarized proton, and the sum of both 
(right). The upper (lower) row gives the results for up (down) quarks. 

are studied as functions of the transverse position and different quark and proton polariza- 
tions taking advantage of the overlap representation of LCWFs that was originally proposed 
in Refs. [HI [12] and successfully applied to GPDs [131 [IS [HI [HI [H] ■ In such an approach 
the lowest-order Fock-space components of LCWFs with three valence quarks are directly 
linked to wave functions derived in constituent quark models. Assuming SU(6) symmetry 
and separating the spin-isospin component $({Aj}, {tj}) from the momentum-space part 
ip{{ki}) of the wave function, the eigenfunctions \I'(r, {Aj}, {tj}) of the light-front Hamilto- 
nian of the nucleon described as a system of three interacting quarks (with helicity Aj and 
isospin Tj, i = 1,2,3) are written as 



*(r,{A,},{r,}) = 2(27r)= 



1 UJiUJ2i^3 



nl/2 



X E ^l{l*(^c/(fci))Dlg:(i?.Kfc2))/^lg;(i?c/(fc3)) $ a,.({A:}, {r,}), (6) 



where r collectively denotes the set of light-cone momentum variables {fc, = {yi,k±i)} 
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of the involved quarks, Mq is the mass of the non-interacting three-quark system, uoi = 
{kf + kl)/^/2, and the matrix D]/^{Rcf{k)) is given by the representation of the Melosh 
rotation Rcf in spin space transforming the canonical spin into the light-front spin. Here the 
momentum space wave function ^/'({fcj}) has the same power- law behaviour with the same 
parameters used in the relativistic quark model of Ref. [TH]. As such the original (instant 
form) wave function only contains S-wave components. 

The various GDPs are then calculated with the above LCWFs as explained in Refs. [T3l 
HU |15] , Fourier transformed according to Eq. ^ and integrated over x to obtain the first 
moments ([S]). 

With such a model only valence quarks are considered. Therefore, the integrals in Eq. ([s]) 
are restricted to < a; < 1. In order to have a first idea of how spin distributions look like 
in a nucleon this limitation is not dramatic because valence quarks are known to dominate 
at large and intermediate x {x > 0.2). 

All the different structures appearing in Eqs. ([T]) and ^ are discussed in the following 
by considering the results shown in Figs. [T] to [7} 

In Figs. [T] and [2] the distorting effect of the dipole terms due to the transverse spin 
distributions on the monopole terms corresponding to spin densities for unpolarized quarks 
in an unpolarized target is shown. If one multiplies the up and down monopole terms by the 
quark charge Cq and sums over quark flavors, one obtains the nucleon parton charge density 
in transverse space. In particular, we find that the central charge distribution for neutron 
target is negative, in agreement with the observations of Ref. [T9] . 

For transversely polarized quarks in an unpolarized proton the dipole contribution intro- 
duces a large distortion perpendicular to both the quark spin and the momentum of the pro- 
ton (Fig. [T]). Evidently, quarks in this situation also have a transverse component of orbital 
angular momentum. This effect has been related [20l [21] to a nonvanishing Boer-Mulders 
function [22] hj- which describes the correlation between intrinsic transverse momentum and 
transverse spin of quarks. Such a distortion reflects the large value of the anomalous tensor 
magnetic moment kt for both flavours. Here, = 3.98 and k,^ = 2.60, to be compared 
with the values /t^ ~ 3.0 and ~ 1.9 of Ref. [9] due to a positive combination + 2H'j'. 
Since kt ~ — ^i", the present results confirm the conjecture that hi is large and negative 
both for up and down quarks [201 EI] ■ 

As also noticed in Refs. [6l [9] the large anomalous magnetic moments k^''^ are responsible 
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Figure 3: The monopole contribution ^SxSx{Ht — A^Ht / 4:M'^) (left) and the quadrupole contri- 
bution ^SxSx{bx — by)H^/M'^ (right) for x-polarized quarks in a nucleon also polarized along x. 
The upper (lower) row gives the results for up (down) quarks. 

for the dipole distortion produced in the case of unpolarized quarks in transversely polarized 
nucleons (Fig. [2|. With the present model, = 1.86 and k*^ = —1.57, to be compared 
with the values = 1.673 and n'^ = —2.033 derived from data. This effect can serve 
as a dynamical explanation of a nonvanishing Sivers function [23] which measures the 
correlation between the intrinsic quark transverse momentum and the transverse nucleon 
spin. The present results, with the opposite shift of up and down quark spin distributions 
imply an opposite sign of for up and down quarks pH [25] as confirmed by the recent 
observation of the HERMES collaboration [26] . 

The results in Figs. [T] and [2] are in qualitative agreement with those obtained in lattice 
calculations [9J, where strongly distorted spin densities for transversely polarized quarks in 
an unpolarized nucleon have been found. One observes that the sideways distortion for down 
quarks is about twice as strong as for up quarks, even if the anomalous magnetic moment 

and the anomalous tensor magnetic moment have about the same magnitude. This is 
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Figure 4: The total spin distribution as a sum of monopole, dipole and quadrupole terms, for 
x-polarized quarks in a proton also polarized along x; left (right) panel for up (down) quarks. 

because the monopole distribution for up quarks is twice as large as for down quarks, and 
therefore adding the dipole contribution results in a larger distortion for down quarks than 
for up quarks. 

In Fig. |3] quarks and proton all transversely polarized along x are considered. Quite 
remarkably, opposite signs of the up and down quadrupole and spin dependent monopole 
terms are found as a consequence of the opposite sign for up and down quarks of the x 
dependence of Ht and Ht predicted by the model [IS]. The quadrupole distribution for up 
quark is more spread than for down quark, but it is much smaller, resulting in an average 
distortion equal to -0.04 to be compared with the value 0.07 for down quark. 

The total spin distribution for quarks and proton transversely polarized along x is shown 
in Fig. |4]as the result of summing for each flavour the two monopole contributions in the left 
panels of Figs. [T] and |3} the two dipole contributions on the middle panels of Figs. [T] and [2] 
and the quadrupole contribution of the right panel in Fig. |3| For up quarks this distribution 
is almost axially symmetrical around a position slightly shifted in the y direction (Fig. |4]). 
This is a consequence of the dominating role of the monopole terms, the compensating action 
of the two dipole terms and the small quadrupole contribution. In contrast the down quark 
spin distribution has a much lower size and shows a strong and symmetric deformation about 
the y axis stretching along the same x direction of the quark and proton polarization. 

The resulting transverse shift of both up and down distributions in Fig. |4] is suggesting 
the presence of an effective transverse quark orbital angular momentum introduced in the 
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Figure 5: The dipole contribution ^SybxE' (left), the total dipole contribution ^[Syb^E' —Sxby{E'rp+ 
2H'rp)/M] (middle) and the quadrupole contribution SxSybxbyH^/M'^ (right) for x-polarized quarks 
in a nucleon transversely polarized in the y direction. The upper (lower) row gives the results for 
up (down) quarks. 

LCWFs by the Melosh rotation required to transform the canonical spin to the light-front 
spin. Due to the shift in the positive y direction transverse quark spin and transverse quark 
orbital angular momentum seem to be aligned along the same x direction of the proton 
polarization. 

In Figs. [5] and |6] results are given for x-polarized quarks in a proton polarized along y. 
The distortion due to the dipole contribution ^SybxE' in Fig. [sjis rotated with respect to 
the case shown in Fig. [TJ but the origin of opposite shift for up and down quarks is always 
the opposite sign of the anomalous magnetic moments /t"''^. The total dipole distortion in 
Fig. [sjis obtained by considering also the second dipole term —^Sxby{E!p + 2H^) /M displayed 
in Fig. [l| The result is quite sizable, while the quadrupole term rather small. Therefore, 
the total resulting distortion of the spin density (Fig. |6]) is due to the dipole terms with a 
small contribution from the quadrupole terms. Correspondingly, the quark orbital angular 
momentum has positive x and y components for up quarks, and positive x and negative y 
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Figure 6: The total spin distribution, as a sum of monopole, dipole and quadrupole terms, for 
x-polarized quarks in a proton transversely polarized in the y direction; left (right) panel for up 
(down) quarks. 

components for down quarks. Here as well as in Fig. |4] the quark orbital angular momentum 
is entirely generated by the Melosh rotations. 

Finally, the case of quark polarization parallel to the proton helicity is considered in 
Fig. [7j Here only monopole terms occur (see Eq. ([T])) and their role was first discussed in 
Ref. [7j. The opposite sign oi \H for up and down quarks is responsible for quite a different 
radial distribution of the axially symmetric spin density. Since in the forward limit the 
GPD H reduces to the helicity distribution gi{x)^ H{x,^ = 0,t = 0) = gi{x), this difference 
ultimately reflects the opposite behaviour of the helicity distributions and the opposite sign 
of the axial- vector coupling constants g^'^ of up and down quarks (see also Ref. p3]). 

In the present analysis only distributions for x > have been considered. In principle, 
also a; < and antiquarks contributions should be considered. However, at the hadronic 
scale of the model contributions from antiquarks are expected to be small in general. For the 
transverse distributions one also knows that gluons do not mix under evolution. Therefore 
their qualitative behaviour presented here in agreement with results from lattice calcula- 
tions [9] can be considered sufficiently indicative. In contrast, sea quarks and gluons can 
affect longitudinal distributions under evolution. In this case, the present results could be 
improved by considering the meson cloud surrounding the bare nucleon at the hadronic scale 
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Figure 7: The monopole contribution t^H (left) and t^H (right), and their sum corresponding to the 
spin distribution for quark polarized in the longitudinal direction, parallel to the proton helicity. 
The upper (lower) row gives the results for up (down) quarks. 

and including its contribution in the evolution [27]. 
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